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Abstract 

We propose the use of Heaviside transform with respect to the quark mass to investigate 
dynamical aspects of QCD. We show that at large momentum transfer the transformed 
propagator of massive quarks behaves softly and thus the dominant effect of explicit chiral 
symmetry breaking disappears through Heaviside transform. This suggests that the mass- 
less approximation would be more convenient to do in the transformed quantity than in 
the original one. As an example of explicit approximation, we estimate the massless value 
of the quark condensate. 



1 Introduction 

In this paper we report the result of the application of Heaviside transform 1 with respect 
to the mass to QCD. Let m be the masses of up and down quarks and the massless value 
of some quantity, Q(m), be of interest. Ordinary approach is to consider the limit that 
lim m ^ Q(m). However, in the framework of perturbation expansion, the limit leads to the 
trivial results and if m is kept non-zero for the approximation of the massless value, the 
effect of explicit chiral symmetry breaking remains and gives, for example, the hard high 
energy behavior of the effective quark mass. In the present paper, we demonstrate that 
the use of Heaviside transform with respect to m solves the above dilemma and provides 
us a new way of approximating the massless dynamics. 

2 Heaviside transform 

Since our approach is based on Heaviside transformation of perturbative series, let us 
state some basic features of the transform. Heaviside transform of Q(m) is given by the 
Bromwich integral, 

~ dm exp(m/m) 

U(m) = / \l\m , 1 

where the vertical straight contour should lie in the right of all the poles and the cut of 
Q(m)/m. Since the Q(m)/m is analytic in the domain, Re(m) > s, £l(rh) is zero when 
m < 0. The Heaviside transform is the inverse of the (second kind of ) Laplace transform 
given by 

f°° - 1 

f2(m) = m j dxexp(— mx)tt(l/x), x — — . (2) 

Since Cl(m) = for rh < 0, the integral (2) reduces to the familiar one in which the 
integration region is [0, oo). However the mathematical manipulation is straightforward 
for the form (2) because Cl(rh) involves Dirac 5 function in usual cases. 
From (2) we find that 

lim Q(m) = lim Cl(rh), (3) 

if the both limits exist. Eq.(3) states that the massless value is an invariant of the Heaviside 
transform. Hence the transformed function can be directly used to approximate Q(0) 
without integration over rh. Here suppose that Q(m) denotes the condensate for the 
quark with its explicit mass m. We note then that, since the massless limit in the both 
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functions cannot be taken at finite orders of perturbation expansion due to the bad infra- 
red behavior of perturbative QCD, we must choose some non-zero m or rh to approximate 
Q(0) or Cl(0), respectively. Then, some difference of the status of approximation may occur 
between the two functions. We will show that, for the approximate calculation of quark 
propagator and the quark condensation, the transformed functions are more convenient 
than the original ones. Though Q is not a physical quantity, there is no objection to use 
it in the approximation of f2(0). 

For the sake of later argument let us show transformations of typical functions. Using 
(1) it is easy to find that 

m -P $ -™L^o{llm) ( P >-1), (4) 



where 9{x) is the step function, 



and 7i represents the Heaviside transformation. By expanding both sides of (4) in powers 
of p, we have 

1 3 9(x), (6) 
\ogm — > -(7 + \ogx)9(x), (7) 

and so on where x = 1/rh. From these results we have the transform of typical functions 
appearing in perturbative expansions. To obtain the formulae, the following result is useful, 

n dCl(l/x) 

milym) — > . (8) 

ox 

For example, using (6) and (8) we find 



m 



k H 



5^\x), (fc = 1,2,3,-..). (9) 



From (7) and (8) we also find 



H 1 

mlogm — > — 9{x) — (7 + \ogx)5(x), (10) 

x 

1 2 / 
m 2 logm — > -^O(x) S(x) — (7 + \ogx)5 (x), (11) 

X 2 X 

2 3 3/ // 
m 3 logm — > ^9(x) -\ — ^5(x) 5 (x) — (7 + \ogx)5 (x). (12) 



The terms containing 8 functions play the role of "counter terms" since they cancel out 
the divergences coming from the Laplace integration of first terms. 

In this paper we confine ourselves with studying the approximate calculation of QCD 
quantities at to — 0. From now on, we therefore omit the 5 function terms of the trans- 
formed functions since the Laplace integration is un-necessary in our scheme that Q(0) will 
be approximated by £l(m) at some fixed non-zero to. 

The result (9) shows that Heaviside transform forces the polynomial of to vanish. This 
is desirable because the polynomial of to should not have nothing to do with the massless 
physics. On the other hand the contributions of the mass-logarithms remains. Since 
the mass-log involves /i, the renormalization parameter, we find that Heaviside transform 
keeps the effect of ultraviolet singularity. This is in accord with the point of view that 
the ultraviolet divergence plays the central role in the non-perturbative effects. Strictly 
speaking, renormalizat ion-group invariant mass, A m , corresponding to to should be used in 
place of to in (1). However for notational simplicity we use to as long as to is proportional to 
A m . We may point out that to obeys the renormalization group equation same as that for 
to. This is because the argument of the exponential in (1) and (2) must be renormalization 
group invariant (This is obvious if one writes transformations (1) and (2) in terms of A m 
and the corresponding conjugate which is necessarily renormalization group invariant). 

In what follows we first study the high energy behavior of the Heaviside function of the 
effective mass and find that it behaves softly for non-zero to at large momentum transfer. 
This fact suggests an advantage of dealing with the Heaviside function when one carry 
out the approximation of the massless case. Next, as an example, we perform a rough 
estimation of the massless value of the quark condensate. Through out this paper we use 
dimensional regularization 2 and Landau gauge. 

3 High energy behavior of the transformed function of the effective quark mass 

The inverse of the quark propagator is written as, 

Sp 1 = p — m — E(to,p) = A(p,m)p — B(p,m) = A(p,m)(p — A4(p, to)), (13) 
and the function M. (p, to) defines the effective quark mass. Heaviside transform of the 
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inverse propagator is given as 

S F 1 = A(p,m)p-B(p,m), (14) 

where 

Sp 1 {p,m) ™ S^fam), A(p,m) % A(p,m), B(p,m) ^ B(p,m). (15) 
The transformed function of the effective mass is defined by 

M( P ,m) = jg4 L (16) 
B{p, m) 

It is easy to show that 

M(p,0) = M(p,0). (17) 

To clarify the difference between A4(p, m) and jti(p,rh), let us discuss the high energy 
behavior of M(p, rh) at the one-loop level. 

At the one-loop level, the self energy of quarks is given by 

-i(g^fC F m(D-l)j^^ £dx[m 2 x-p 2 x(l-x)Y\ (18) 

where Cf = (N 2 — 1)/2N C and D = 4 — 2e. Then at large — p 2 the inverse of quark 
propagator behaves as 

S" 1 = p-M(p,m), (19) 
MM = m l-_(log-^ --) +-— (log ^ + l)+0(^), 



where we have subtracted the ultraviolet divergence according to MS scheme 3 which we 
use through out this paper. Heaviside transform of (19) will be performed by using the 
following results coming from (6), (9) and (12), 

1^1, m ™ 0, m 3 (logm) ^ -2m 3 . (20) 

We note from (20) that the so-called hard piece of order m in (19) disappears and the soft 
term of order m 3 /p 2 survives after the transform. The reason that the soft term remains is 
that it involves log(m). The fact that the mass-log exists at the order m 3 /p 2 can be found 
by differentiating (18) with respect to m. After the differentiation up to three times, one 
finds the infra-red singularity when m — > 0. 



From (19) and (20) we have 

S- 1 = p-M(p,m), (21) 
M(p,m) = +0(1 p*), 

7T p A 

and find that Ai behaves softly []. We note that, for non zero rh, Ai(p,rh) simulates the 
soft behavior of Ai(p, 0) with spontaneous symmetry breaking 5 . Thus, the dominant effect 
of explicit chiral symmetry breaking has been washed out through Heaviside transform. 
From the second equation in (20) the hard piece vanishes to all orders if it is simply 
proportional to m. Actually it is pointed out in Ref.6 that, in pole subtraction schemes, 
the coefficient functions of operator product expansion 7 are analytic in all masses involved. 
Then, since the hard piece belongs to the unit operator and the unit operator cannot 
involve any mass-log, the whole part of the order m is just proportional to m. Thus, we 
conclude that the leading effect of the explicit mass disappears in the Heaviside function 
of the effective mass to all orders. This reveals an advantage of Xi (p, rh) in approximating 
M.(p, 0), since keeping rh non- zero is compatible with the chiral symmetry. On the other 
hand, for Ai(p, m), we must set m = to force the hard piece vanish but this makes the 
perturbative M. trivial. 

To carry out the massless approximation, we must take higher order contributions into 
consideration and approximate M.{p, 0) by setting rh as small as possible in M.(p, rh) |. A 
typical example is given in the next section. 

4 Calculation of the quark condensate 

To prove non-perturbative effects, the perturbative series would not contain enough 
information. However, it is not clear whether it is useless in the approximate calculation 
or not. Actually we have found that the perturbative series is effective in the approxi- 
mate calculation of the dynamical mass in the Gross- Neveu model 8 . Since the Heaviside 

transform removes the dominant effect of the explicit mass which is irrelevant to massless 

tThe soft behavior was also found by the direct calculation in the generalized Hartree-Fock approach 4 . 
Also, we point out that, though all sub-leading terms in l/p 2 are lost at the one-loop level, higher order 
contributions would involve log to and give non-trivial results under TL transformation. 

t Higher order contributions would produce powers of log to at the rh?/p 2 order term in Ai and the 
resulting series of to 3 times the series of \ogrh/n corresponds to the transformed quark condensate. 
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dynamics, we expect also in QCD that the transformed quantity is more convenient in the 
approximate calculation than the original one. 

In this section we try to estimate the massless value of the condensate. For the purpose 
we directly calculate the condensate for the massive case to two-loops and transform it to 
evaluate the massless value, assuming the existence of the non-zero condensate. 

As is well known, the naive product, J2a,o q A ' a (fyq A (fy = 59> * s singular (A and a denote 
color and Lorentz indices, respectively) and we need to regularize it. In this paper, we 
define the regular product according to MS scheme 9 . We denote thus regularized product 
as [qq]. In this definition the two- loop calculation of the condensate is given by 10 
., .. N c m 3 u- 2e (^ , m 2 3C F g 2 / 2 m 2 5 , m 2 \\ 

» = -sH 1 - log u + ^( log i? - 3 log 7? + constant ) )• (22) 



where 

r i r, N c m 3 a- 2e (\ C F g 2 , \ 2 ix\ . . 

[*] = z *~t?-U + -^-( 3 7-7)> (23) 

and 

3C F g 2 (fi)i 11 . . 

Note that we need non-multiplicative renormalization to make the condensate finite (see 
(23)). As a result, the transformation property of the product under the renormalization 
group changes from that of Zqq and {[qq]) no longer satisfies /id{[qq]) / dfj, = — 7 m {[qq]) 
where 7 m m = udm/dfi = (— ^§^a + 0(a 2 )^m. 

Our task is first making Heaviside transform of (22) and then improving the result by 
renormalization group. First, we point out that, since the non-multiplicative pole terms 
are the cube of m, they disappear after Heaviside transform to give, 

([¥]) = {Zqq), (25) 

where the hat denotes the Heaviside transform and we used (see (9)) 

m — > (k = positive integer). (26) 

Since in the pole subtraction scheme the non-multiplicative piece is proportional to m 3 
and free from the mass-log, eq.(25) holds to all orders. Thus, as in the exact massless 
case, we could use Zqq as the local product. As well as the recovery of the soft behavior 



in the transformed effective mass, this result shows that the effect of the explicit mass is 
reduced and it would be better to use the transformed quantity to simulate the massless 
case. Now, from (9), (12) and 



m 3 log 2 (m) — 2(m/n) 3 (3 + 2(k>g(m/n) — 7) 



we find 



H /f r- n N c m A 



« - (Mi = 



7T Z 



1 



3Cpa(fj J ) 



IT 



, m 2 13\ 



(27) 



We find by the direct operation of fi(d/dfi) that the condensate satisfies the renormalization 
group equation, 

= -7m([gg]), 7™ = — + 0(a 2 ), (29) 

up to the lowest order. Thus, as implied by (25), the transformed condensate shows the 
correct property under the renormalization group. The solution of (29) is given by 

9C F 



([??]) = X{a,rh, {M'},/i)| M=Mo a(^) , X 



a A A 



(30) 



\\C G -2n F 

where the subscript means the value at /i — /i and {M'} denotes the set of explicit 
masses of other quarks. 

Let us improve the large m behavior of x which can be settled by perturbation expan- 
sion. The behavior is improved by adjusting /i in accordance with m . We impose /io to 
satisfy 

log— -7 + — = t, (31) 
/io iz 

where t is a fixed constant which value is yet un-specified. This condition fixes /i as the 
function of A m (the renormalization group invariant mass corresponding to m), A (finite 



QCD scale in MS scheme), and t. Further the A m dependence of mo changes from the 
simple proportional one. Actually, from the renormalization group equation and (31), we 
find the implicit equation for m , 

-7+13/12-t\ ~ A 



m = A m a(/i ) =Am[Po log 



m e 



A 



(32) 



at this order of expansion. Since /io > A to keep the coupling a positive, rho must be 
larger than Ae +7 ~ 13 / 12+ *(= to*) from (31). When A m goes to zero, m — > rh* and a — > +00 



from (31) and (32). Thus, rh* represents the limitation of the obtained perturbative result. 
Although the correspondence between A m and m has been involved, rh is a monotonic 
function of A m and most manipulation can be carried out in terms of m . We note from 
(31) that a depends on m as governed by rhodao/dm = (3(a ) = — Poa^ + 0(al) where 
(3 = (IICg — 2ni?)/67r. The asymptotic freedom in QCD 11 enters into the condensate 
through ao(m ,t). Taking (31) into account, we thus have 



N c m n 

and arrive at the improved result, 



3 r 



6Cj?a!o(mo> t) 



TT 



a (rh ,t) A , (33) 



N r fhl r i _ 6C F ao(m ,t) ^ 

TT 



-A _ ^C'- Q 
TT 2 



a (m ,t) A ■ a(fi) A . (34) 



Let us turn to the rough estimation of the condensate, ([qq]) | m =o- Recall that one cannot 
take the A m — > limit as we stated in the previous paragraph. Then what we can do is to 
fix A m or equally rho in terms of A under the guiding principle that one should minimize the 
effect of the explicit mass as possible as one can. Here we note that linim^o m? ( [qq] ) = 
for any posotive integer j and therefore from (3) that lim ri j^ H[m J ([??])] = (H denotes 
the operation of Heaviside transformation). Let us impose the case of j = 1. Further, 
since we have another parameter t which we can choose freely, we demand the condition 
of j = 2. Now from (8), the conditions of rh Q minimizing the effect of the explicit mass are 

d(l/m ) ' d(l/m ) 2 ' { } 

which lead to the same conditions for x- Eq.(35) gives two non-trivial solutions. However, 
one of those corresponds to the large coupling, a ~ 8, and we discard this since it is 
too large to rely upon within the two-loop level. The other solution gives a ~ 0.8 for 
various flavors and is valuable to be discussed further. That value of a indicates that the 
renormalization scale /i is a few of A. Assuming the decoupling of heavy quarks 12 , we 
therefore work with the choice that n_p = 3. Then we have a ~ 0.79 and t ~ 0.91. These 
values give m ~ 3.6A 3 , x ~ — IO.8A3 and A m ~ 4A 3 where (31) and the one-loop relation, 
A 3 = /iexp(— l/(3 a), was used (A 3 denotes the QCD scale effective at three flavors). Thus 
we arrive at 

([¥]) ~ -10.8A^- A (/x). (36) 



Further numerical computation of the condensate can be done as follows: From the 
experimental data, a(mz) ~ 0.117 ± 0.005 13 , and the definition of A at the lowest order, 
we have A 5 ~ 83 + 27 — 23MeU for n F = 5. With the help of the relation, A 3 ~ 
A 5 (m c m b /Al) 2 / 27 (m b = A.3GeV,m c = 1.3GeU) 14 , it is converted to the scale value for 
tif = 3, giving A 3 ~ 135 + 28 — 3lMeV. By setting fj, = lGeV and n F = 3 in (36), we 
thus have 

(\m\)\n=iGeV ~ -(351 + 88 - 77 MeVf. (37) 

5 Discussion 

The phenomenological value, (uu) ~ (dd) ~ — (250 MeV) 3 15 is slightly below the range 
shown in (37). We think that the result is rather good, taking into account that the result 
(37) is just the lowest order one which comes from a rough ansatz, (34). Of course the 
higher order calculation is needed for the serious study of our approach. In particular the 
three loop calculation is important since at this order the non-linear gluon vertex and the 
loops of other quarks becomes explicitly active. 
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